clear all;
clc;
folder='sigma2_29Apr08\';
folder2='sigma2_30Apr08\';
%
tmax=1000;
T=1:tmax;

sig2_0_025 = dlmread(strcat(folder,'sig2_0_025.dat'));
sig2_0_025_up = dlmread(strcat(folder,'sig2_up_0_025.dat'));
sig2_0_025_down = dlmread(strcat(folder,'sig2_down_0_025.dat'));


sig2_0_05 = dlmread(strcat(folder,'sig2_0_05.dat'));
sig2_0_05_up = dlmread(strcat(folder,'sig2_up_0_05.dat'));
sig2_0_05_down = dlmread(strcat(folder,'sig2_down_0_05.dat'));


sig2_0_1 = dlmread(strcat(folder,'sig2_0_1.dat'));
sig2_0_1_up = dlmread(strcat(folder,'sig2_up_0_1.dat'));
sig2_0_1_down = dlmread(strcat(folder,'sig2_down_0_1.dat'));


sig2_0_2 = dlmread(strcat(folder,'sig2_0_2.dat'));
sig2_0_2_up = dlmread(strcat(folder,'sig2_up_0_2.dat'));
sig2_0_2_down = dlmread(strcat(folder,'sig2_down_0_2.dat'));

sig2_0_4 = dlmread(strcat(folder2,'sig2_0_4.dat'));
sig2_0_4_up = dlmread(strcat(folder2,'sig2_up_0_4.dat'));
sig2_0_4_down = dlmread(strcat(folder2,'sig2_down_0_4.dat'));

sig2_0_8 = dlmread(strcat(folder2,'sig2_0_8.dat'));
sig2_0_8_up = dlmread(strcat(folder2,'sig2_up_0_8.dat'));
sig2_0_8_down = dlmread(strcat(folder2,'sig2_down_0_8.dat'));

figure(1)
plot(T,sig2_0_025,T,sig2_0_05,T,sig2_0_1,T,sig2_0_2,T,sig2_0_4,T,sig2_0_8);
grid on;
xlabel('T:1000');
ylabel('{\sigma^2}_\epsilon(t)');
title('{\sigma^2}_\epsilon(t) vs T for \epsilon = 0.025, 0.05, 0.1 and 0.4,0.8 for 100 iterations');


% figure(2)
% plot(T,sig2_0_025,T,sig2_0_025_up,T,sig2_0_025_down);
% grid on;
% xlabel('T:1000');
% ylabel('{\sigma^2}_\epsilon(t) for charge, up and down');
% title('{\sigma^2}_\epsilon(t) charge, up, down vs T for \epsilon = 0.025 for 200 iterations');

% figure(3)
% plot(T,sig2_0_05,T,sig2_0_05_up,T,sig2_0_05_down);
% grid on;
% xlabel('T:1000');
% ylabel('{\sigma^2}_\epsilon(t) for charge, up and down');
% title('{\sigma^2}_\epsilon(t) charge, up, down vs T for \epsilon = 0.05 for 200 iterations');
% 
% 
% 
% figure(4)
% plot(T,sig2_0_1,T,sig2_0_1_up,T,sig2_0_1_down);
% grid on;
% xlabel('T:1000');
% ylabel('{\sigma^2}_\epsilon(t) for charge, up and down');
% title('{\sigma^2}_\epsilon(t) charge, up, down vs T for \epsilon = 0.1 for 1000 iterations');
% 
% 
% figure(5)
% plot(T,sig2_0_2,T,sig2_0_2_up,T,sig2_0_2_down);
% grid on;
% xlabel('T:1000');
% ylabel('{\sigma^2}_\epsilon(t) for charge, up and down');
% title('{\sigma^2}_\epsilon(t) charge, up, down vs T for \epsilon = 0.1 for 1000 iterations');

figure(6)
plot(T,sig2_0_4,T,sig2_0_4_up,T,sig2_0_4_down);
grid on;
xlabel('T:1000');
ylabel('{\sigma^2}_\epsilon(t) for charge, up and down');
title('{\sigma^2}_\epsilon(t) charge, up, down vs T for \epsilon = 0.4 for 100 iterations');


figure(7)
plot(T,sig2_0_8,T,sig2_0_8_up,T,sig2_0_8_down);
grid on;
xlabel('T:1000');
ylabel('{\sigma^2}_\epsilon(t) for charge, up and down');
title('{\sigma^2}_\epsilon(t) charge, up, down vs T for \epsilon = 0.8 for 100 iterations');

% figure(2)
% plot(T,sig2_0_025);
% 
% figure(3)
% plot(T,sig2_0_05);
% 
% figure(4)
% plot(T,sig2_0_1);
% 
% figure(5)
% plot(T,sig2_0_2);
